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Abstract. Early results from the INTEGRAL Core Program, for a sample of eight persistently bright neutron star low mass X- 
ray binaries in the energy range from 5 keV to 200 keV are presented. It is shown that INTEGRAL efficiently detects sources and 
that spectra may be obtained up to several hundreds of keV by combining data from three of the four INTEGRAL instruments: 
JEM-X, IBIS and SPI. For the source GX 17+2 it is shown that the spectrum extends well above 100 keV with a flattening of 
the spectrum above 30 keV. This might suggest a non-thermal comptonisation emission, but uncertainties in the current data 
reduction and background determination do not allow firm conclusions to be drawn. 
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1. Introduction 

Since its launch in October 2002, the International Gamma- 
Ray Astrophysics Laboratory, INTEGRAL, has been providing 
a large amount of interesting data. The combination of the two 
wide field of view (FoV) instruments, the imager IBIS (15 keV 
- lOMeV, 29° x 29° partially coded FoV, Ubertini et al. 2003) 
and the spectrometer SPI (20keV - 8 MeV, 35° x 35° partially 
coded hexagonal FoV, Vedrenne et al. 2003) coaligned with the 
JEM-X (Lund et al. 2003) and OMC (Mas-Hesse et al. 2003) 
monitors, allows large areas of the sky to be observed and mon- 
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itored in one pointing in a wide frequency range from the opti- 
cal to the 7-ray domain. Such a capability is fully exploited dur- 
ing the INTEGRAL Core Program (a series of successive scans 
of the Galactic Plane (GPS; Winkler et al. 2003) and Galactic 
Centre (GCDE)) which is regularly producing large amounts of 
data, in particular on persistently bright sources. 

The aim of this paper is to report preliminary results from 
early measurements on eight persistent bright Low Mass X-ray 
binaries (LMXRB) hosting a neutron star. The sample has been 
selected from a larger set observed during the Core Program 
Observation scans on the Galactic Plane executed so far. The 
sources are listed in Tab. [2 Given the type of sources involved 
(about hundred of mCrabs in the 2-10keV band, Liu et al. 
2001) and the pointing exposures of about 2000 s, for sin- 
gle pointing spectral extraction, JEM-X for soft photons (5- 
20 keV) and the low energy IBIS detector, ISGRI (Lebrun et al. 
2003) for harder photons (20-200 keV) were chosen. PICsIT, 
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the hard photon IBIS detector (Di Cocco et al. 2003), has its 
peak sensitivity above 200 keV while the LMXRB are consid- 
ered to display a rather steep spectrum and consequently have 
fluxes below PICsIT detectability. The spectrometer SPI has 
been used to extract the hard energy spectra (20-200 ke V) av- 
eraged on longer time scales. 

The combination of JEM-X, IBIS and SPI data provides a 
complete soft to hard energy coverage, allowing a regular mon- 
itoring of source behaviour. Special attention is given to the 
hard emission (> 50 keV) monitoring for which INTEGRAL 
will give unprecedented continuous coverage. 

Sec. 2 of this paper gives an overview of the state of 
the art studies on LMXRBs containing a weakly magnetized 
neutron star and a short description of the INTEGRAL Core 
Program selected sample. Sec. 3 contains details of the rel- 
evant INTEGRAL observations and data reduction methods 
used. Preliminary results are summarised in the last Section. 



2. Bright LMXRBs in the INTEGRAL Core 
Program 

2. 1. LMXRBs with weakly magnetised neutron stars 

LMXRBs hosting a weakly magnetised neutron star can be 
broadly classified into two classes (van der Klis, 1995): high 
luminosity/Z sources and Atoll sources covering a much wider 
range in luminosity. Z sources describe an approximate "Z"- 
shape (horizontal, normal, and flaring branch) in the colour- 
colour (CC) and X-ray hardness intensity diagrams while Atoll 
sources are characterised by an upwardly curved branch. Two 
recent studies (Muno et al. 2002 and Gierlinski & Done 2002) 
suggest that the clear Z/Atoll distinction on the CC diagram 
is an artifact due to incomplete sampling: Atoll sources, if ob- 
served long enough, do exhibit a Z shape in the CC as well. 
Many differences, however, remain: Atoll sources have weaker 
magnetic fields (about 10 6 to 10 7 G versus 10 8 -10 9 G of Z 
sources), are generally fainter (0.01-0.3LEdd versus ~ LEddX 
can exhibit harder spectra, trace out the Z shape on longer time 
scales than typical Z-sources and have a different correlated 
timing behaviour along with the position on the Z. Thus the 
distinction, at least from a practical point of view, still makes 
sense. 

Recent broad band studies, mainly with BeppoSAX, have 
shown that many Z sources display a variable hard power-law 
shaped component, dominating their spectra at > 30 ke V (Di 
Salvo & Stella l2002l D' Amico et al. l200"Tl Di Salvo et al.|2000 
and references therein). This power-law has been explained as 
non-thermal Comptonisation. As previously said, lower lumi- 
nosity systems (Atolls) can display much harder spectra which 
can be well described with a simple power-law with photon 
indices of about 1.5-2.5. Hard X-ray components extending 
up to a few hundred keV have been seen in about 20 neutron 
star LMXRB of the Atoll class. These sources usually dis- 
play an exponential cut-off between many tens and a few hun- 
dred keV. This component is interpreted as unsaturated ther- 
mal Comptonisation and is known as the "hard state" of Atoll 
sources. 
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Fig. 1. Exposure map for a total of 6 months of GPS and GCDE data. 
The spatial distribution of the 8 sources is visible (red crosses). 



The long term X-ray variability of those sources has been 
extensively studied in the 2-12 keV band with the RXTE All 
Sky Monitor, as well as during RXTE dedicated pointings till 
about 40 keV (van der Klis l2000l Swank |200B . On the other 
side, BeppoSAX pointings have shown the presence of the hard 
tails mentioned above showing that Neutron Star systems, as 
well as Black Hole ones, are capable of producing such hard 
photons. 

The combination of regular monitoring in the hard X- 
rays and 7-rays has not been done before and this is where 
INTEGRAL will give a major contribution to understanding the 
behaviour of bright Low Mass X-ray binaries from 5 keV up to 
about 200 keV. 

2.2. The INTEGRAL Core Program sample 

Thirtyfive per cent of INTEGRAL observing time, the Core 
Program, is time reserved for the institutes that developed and 
delivered the instruments, for the INTEGRAL Science Data 
Centre (ISDC; Courvoisier et al. 2003) and for the Russian 
scientists in return for providing the Proton rocket which put 
INTEGRAL in orbit. 

Fig. 1 shows the exposure map for a total of six months 
of this program. As can be seen, the centre of our Galaxy has 
been heavily covered, with the total observing time decreasing 
as we move away from it. From the overall list of about 70 per- 
sistent and transient LMXRBs, the 8 sources detected primar- 
ily by IBIS/ISGRI and then by JEM-X in at least 20 pointings, 
have been selected. At the time of writing OMC data on these 
sources is still quite sparse and, therefore, will not be included 
in this paper. The red crosses in Fig. 1 show the position of 
those sources in the Galaxy while in Table 1 the complete list 
is presented. 

The exposure per source depends on its position relative to 
the scan path and will differ from instrument to instrument due 
to their different FoV. 
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Name 

Sco X-l 

1822-371 

GX3+1 

GX9+9 

GX9+1 

GX5-1 

GX 17+2 

CygX-2 



1 

359.09 
356.85 
2.29 
8.51 
9.08 
5.08 
16.43 
87.33 



23.78 
-11.29 
0.79 
9.04 
1.15 
-1.02 
1.28 
-11.32 



Type 

Z 
ADC 
AB 

A 

A 

Z 

ZB 
ZB 



/rxte 
11420± 1875 
25 ± 29 
390 ± 64 
265 ± 44 
495 ± 75 
913 ± 137 
603 ± 104 
440 ± 80 



Table 1. Bright persistent LMXRBs regularly monitored 
by INTEGRAL. Type: A=Atoll, B=bursting, Z=Z-source, 
ADC=Accretion Disc Corona; Flux: average fluxes (in mCrab) 
as observed by RXTE/ASM (1.5 - 12keV) and by SPI (20 - 
40keV; about 1 Msec overall exposure). See text. 




3. Data reduction and analysis 

A large fraction of the Galactic Centre Deep Exposure (GCDE) 
has already been completed. One scan of the GPS is performed 
every 12 days on average. 

We have analysed most of the Core Program data currently 
available with the standard INTEGRAL Data Analysis System 
(IDAS). The following 3 subsections provide the instrument 
specific analysis details. 

3.1. SPI analysis 

The analysis of the SPI data is based on GCDE observations 
taken from revolution 47 up to 62, i.e. between March 3rd and 
April 19th, 2003. 544 dithering pointings used in the analysis 
combine a total exposure time of 958 ksec. As the instrumen- 
tal resolution of SPI is about 2.5°, source confusion can affect 
the results. The ISGRI data have been used as a reference to 
check for sources which might interfere in the SPI data. GX 5- 
1 is separable by 40 arcmin from the black hole candidate GRS 
1758-258. In this case both sources will influence the results of 
source extraction of each other and fluxes and spectra can only 
be taken as tentative. In the vicinity of GX 17+2 two sources 
can be detected, 4U1812-12 and AX J1820.5-1434, at 2.0° and 
1.2° angular separation, respectively. In these cases the SPI 
Iterative Removal Of Sources (SPIROS) program (Skinner & 
Connell, 2003) is able to disentangle the sources fairly well as 
the closer one is rather faint, though minor effects cannot be 
excluded. SAX J 1805 .5-2031 is located 1.0° away from GX 
9+1. 

As both sources appear to have similar fluxes in the studied 
energy range, the spectrum of GX 9+1 derived from the SPI 
data might therefore be affected. In the case of GX 3+1, two 
faint sources, SLX 1735-269 and SAX J1747.0-2858, are at 
2.2° and 2.3° angular separation, respectively. As both sources 
are weak compared to GX 3+1, significant effects on the spec- 
tral extraction of GX 3+1 are unlikely. Fig. |2] shows an im- 
age extracted from the SPI data in the 40-100 keV region. In 
the very dense region of the Galactic Centre, the extraction of 
sources seems to fail, while for more isolated sources, like e.g. 
GS 1826-24 and GX 13+1, the results are consistent with the 
ISGRI data (see Fig. For spectral extraction twenty loga- 
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Fig. 2. SPI image of the Galactic Centre in the 40 - 100 keV 
band. The sources in white belong to our monitoring program. 
Other sources are marked in yellow for orientation. 



rithmic energy bins in the 20-200 keV energy range have been 
applied to the data (Fig. 0). The instrumental response func- 
tion used for the analysis, has been derived from on-ground- 
calibration and corrected after the in orbit Crab calibration ob- 
servation. Source fluxes in the 20^10 keV band have been com- 
puted by comparison with results from Crab observation and 
are therefore given directly in Crab units (see Tab. 0. Fluxes 
in the 1.5-12keV energy band have been extracted from the 
RXTE/ASM data base. The fluxes have been averaged over the 
same time period as the SPI data in order to have comparable 
results. In the case of 4U 1822-371 only the measurements for 
which the flux added to the la error is larger than 0.0 have been 
taken into account, so that only reliable data are recognized. 

3.2. ISGRI 

The analysis of ISGRI data is based on GCDE and GPS data 
from revolution 30 to 64 i.e. January 1 1th to April 22nd, 2003. 
One thousand pointings (for a total of about 2 Msec exposure 1 ) 
have been analysed separately and then combined in the mo- 
saicked image shown in Fig. [3] a zoom in the Center of the 
Galaxy in the 20^10 keV and 40-60keV bands. 

Given the ISGRI sensitivity (5a detection in the 20^-0 keV 
band for a 20 mCrab source in one pointing of 2200 sec; 
Rodriguez et al. 2003, possible systematic errors are not taken 
into consideration) and the brightness of the sources of our 
sample, we can extract a source spectrum using (unlike for SPI) 
only data from one pointing. The ISGRI spectrum for a single 
pointing from revolution 54 for GX 17+2 is shown in Fig. [4] 



1 This is the total exposure time. The single source (point) exposure 
is much less as can be seen in Fig.Q 
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Fig. 4. GX 17+2 combined JEM-X (black), ISGRI (green) and 
SPI (red) count spectra. The SPI spectrum is based on all 
GCDE data from March to April 2003, while JEM-X and 
ISGRI spectra are extracted from one pointing only (about 2.2 
ksec). Individual normalisation has been applied to the three in- 
struments in order to compensate for uncertainties in the cross- 
calibration. 



3.3. JEM-X 

The set of data taken from the same position for the same ex- 
posure time used for ISGRI has also been analysed for JEM- 
X using the standard software. Only JEM-X2 data have been 
available during the observations performed so far. The statis- 
tics for a bright LMXRB like GX 17+2 is of course even better 
than in ISGRI (1 pointing from revolution 54 extraction) and 
the result is shown in Fig.|4] 

Light-curves have also been extracted for each source and 
for each pointing and a 1 day lightcurve for GX 3+ 1 is shown as 
an example in Fig.|5]for two different energy bands (5-12 keV 
and 12-20 keV). For comparison the Crab pulsar shows a count 
rate in JEM-X of 70 — 75 counts s -1 and 20 — 25 counts s _1 in 




Fig. 5. JEM-X light-curve for GX 3+1 in two energy bands. 
Each point refers to a 100 second time bin. IJD is the fractional 
number of days since January 1, 2000 (0 UT). See text for more 
details. 



the 5-12 keV and 12-20 keV energy band, respectively. Gaps in 
the data mean that for those periods the source was out of JEM- 
X FoV, the narrowest of INTEGRAL high energy instruments, 
thus there is no detection. 



4. Results 

Since its launch, INTEGRAL has been producing a huge 
amount of interesting data. Quite naturally at this early stage of 
the mission, however, it is not possible to exploit them at best 
mainly because instrument response functions are still under 
development and calibrations on-going. More time is required 
before one can present long and fast time scale variabilities in 
different energy bands as well as detailed spectral behaviour of 
the sources. Nonetheless some main conclusions can be already 
drawn. First of all, the combined spectrum in Fig.|4]shows the 
INTEGRAL mission wide energy band allows the emission of 
the binary system to be displayed all the way to the hard tail. 
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There seems to be no evidence for a cut-off below 100 keV. In 
this case the combined spectrum can be represented by a simple 
model, including a blackbody plus a comptonisation compo- 
nent. The statistical quality of the fit does not justify adding an 
additional power law component. However, the present model 
flux is much below the measured one in the channel centered 
on 100 keV. This is a possible indication of the presence of a 
high energy tail on top of the blackbody Comptonised by ther- 
mal electrons. Physically, this can be realized if the electrons 
have a hybrid distribution, with a Maxwellian and a high en- 
ergy tail. This appears to be the case in a number of black hole 
binaries (e.g., Gierlinski et al. 1999; Zdziarski et al. 2001). 

Furthermore, based on one pointing only, JEM-X and 
ISGRI give a significant spectrum up to 100 keV for a source of 
~ 60 mCrab. This will allow the study of the spectral time evo- 
lution on a pointing by pointing basis, i.e. monitoring the spec- 
tral slope and cut-off on an hourly basis. In this energy region 
the SPI sensitivity is much lower when compared to ISGRI. 
The spectrograph can then be used to derive high energy spec- 
tra based on longer observations. Though the spectrum is in 
this case an average over the different stages of the LMXRB 
cycle, it can reveal high energy tails, as seen in Fig. |4]up to 
~ 150 keV. This capability is also seen in the high energy (40 
- 100 keV) image as shown in Fig.|2] 

In addition, note that the ISGRI FoV is about the size of the 
image in Fig. [5] so the amount of information one can obtain 
from one ISGRI pointing is evident. Besides, its excellent an- 
gular resolution (12' FWHM) allows the emission from close- 
by sources to be separated. See for instance GX 5-1 and the 40 
arcmin distant Black Hole Candidate GRS 1758-258, for which 
possible SPI source confusion has been pointed out. 

Another result of this study is the comparison of the SPI 
derived fluxes to the RXTE/ASM fluxes as shown in Tab. [2 As 
both flux values have been extracted over a long time period 
of 1 .5 months, they represent the mean status of the LMXRBs, 
i.e. averaged over all positions in the colour-colour diagram. 
While the Z and Atoll sources of the sample show a steep spec- 
tral slope between the 1.5-12keV and the 20^40 keV band, 
the 0.59 second pulsar 4U 1822-371 (Jonker & van der Klis 
2001 1 exhibits a comparably flat spectrum with a photon index 
of T ~ 2 between the RXTE and SPI data. This hard spectrum 
has been reported before (e.g. Parmar et al. 2000 1, but further 
investigations are necessary in order to understand the nature 
of this illusive source. 

Finally, though the X-ray monitor JEM-X covers a smaller 
area than the two main instruments, Fig.|5]shows that monitor- 
ing of the X-ray flux on short (several minutes) and long (days 
and months) time scale is possible. This will allow to monitor 
outburst and the dipping behaviour of 4U 1822-371. 
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